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Abstract 

This paper is concerned with the pure-state A-representability problem for systems under 
a magnetic field. Necessary and sufficient conditions are given for a spin-density 2x2 matrix 
R to be representable by a Slater determinant. We also provide sufficient conditions on the 
paramagnetic current j for the pair (i?,j) to be Slater-representable in the case where the 
number of electrons N is greater than 12. The case A < 12 is left open. 


1 Introduction 

The density-functional theory (DFT), first developed by Hohenberg and Kohn [T], then further 
developed and formalized mathematically by Levy [5], Valone [3] and Lieb [3], states that the 
ground state energy and density of a non-magnetic electronic system can be obtained by mini¬ 
mizing some functional of the density only, over the set of all admissible densities. Characterizing 
this set is called the N-representability problem. More precisely, as the so-called constrained 
search method leading to DFT can be performed either with A-electron wave functions nil], 
or with A-body density matrices mill, the A-representability problem can be recast as follows: 
What is the set of electronic densities that come from an admissible N-electron wave function? 
(pure-state A-representability) and What is the set of electronic densities that come from an ad¬ 
missible N-electron density matrix? (mixed-state A-representability) This question was answered 
by Gilbert [S], Harriman [3| and Lieb [I] (see also Remark [T|). 

For a system subjected to a magnetic field, the energy of the ground state can be obtained by 
a minimization over the set of pairs (R, j), where R denotes the 2x2 spin-density matrix [7] (from 
which we recover the standard electronic density p and the spin angular momentum density m) 
and j the paramagnetic current [8|. This has lead to several density-based theories, that come from 
several different approximations. In spin-density-functional theory (SDFT), one is only interested 
in spin effects, hence the paramagnetic term is neglected. The SDFT energy functional of the 
system therefore only depends on the spin-density R. The A-representability problem in SDFT 
are therefore: What is the set of spin-densities that come from an admissible N-electron wave 
function? (pure-state representability) and What is the set of spin-densities that come from an 
admissible N-hody density matrices? (mixed-state representability). This question was left open 
in the pioneering work by von Barth and Hedin [^, and was answered recently in the mixed-case 
setting [Tj. In parallel, in current-density-functional theory (CDFT), one is only interested in mag¬ 
netic orbital effects, and spin effects are neglected m- In this case, the CDFT energy functional 
of the system only depends on p and j, and we need a characterization of the set of pure-state and 
mixed-state A-representable pairs (p,j). Such a characterization was given recently by Hellgren, 
Kvaal and Helgaker in the mixed-state setting El, and by Lieb and Schrader in the pure-state 
setting, when the number of electrons is greater than 4 m- In the latter article, the authors 
rely on the so-called Lazarev-Lieb orthogonalization process [T2| (see also Lemma [S|) in order to 
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orthogonalize the Slater orbitals. 

The purpose of this article is to give an answer to the A^-representability problem in the 
current-spin-density-functional theory (CSDFT): What is the set of pairs that come from 

an admissible N-electron wave-function? (pure-state) and What is the set of pairs that 

come from an admissible N-body density-matrix? (mixed-state). We will answer the question in 
the mixed-state setting for all iV e N*, and in the pure-state setting when > 12 by combining the 
results in [7] and in [12]. In the process, we will answer the Wrepresentability problem for SDFT 
for all N in the pure-state setting. The proof relies on the Lazarev-Lieb orthogonalization 
process. In particular, our method does not give an upper-bound for the kinetic energy of the 
wave-function in terms of the previous quantities (we refer to [Idl I14| for more details). We leave 
open the case < 12 for pure-state CSDFT representability. 

The article is organized as follows. In Section [51 we recall briefly what are the sets of interest. 
We present our main results in Section jS] the proofs of which are given in Section |3| 


2 The different Slater-state, pure-state and mixed-state sets 

We recall in this section the definition of Slater-states, pure-states and mixed-states. We denote 
by C'“°(M^), ... the spaces of real-valued L^, (7““, ... functions on R.^, and 

by LP(R^C‘^), C'“(R3^C'^), ... the spaces of C^-valued LP, H\ C°° functions on 

K^. We will also make the identification L^’(IR^,C‘^) s (LPfR^,<C)Y (and so on). The one-electron 
state space is 

l2(m3^C') e 1$ = ll‘i>IU= == Xs ’ 

endowed with the natural scalar product (<I>i|$ 2 ) ■= /rs +'/’i</' 2 )- "^^^e Hilbert space for N- 

electrons is the fermionic space A^i T^(M^, C^) which is the set of wave-functions e 
satisfying the Pauli-principle: for all permutations p of {1,..., A^}, 

■ • ■ 7 ^p{N) ^ ^p{N)') ~ ■ • ■ 7 7 ^A^") 7 

where e{p) denotes the parity of the permutation p, e position of the fc-th electron, and 

Sk 6 {t 7 -l} its spin. The set of admissible wave-functions, also called the set of pure-states, is the 
set of normalized wave-function with finite kinetic energy 

:= jvl/ 6 A l 2(R3,C2), lIVvI-llA < 00 , ||VI-||A(R3«) = l| 

where V is the gradient with respect to the 3N position variables. A special case of wave-functions 
is given by Slater determinants: let $i, $ 2 ,..., be a set of orthonormal functions in L^(R^, C^), 
the Slater determinant generated by ($ 1 ,..., $jv) is (we denote by x*, := (r^, Sk) the fc-th spatial- 
spin component) 

^[$i,...,$Ar](xi,...,xjv) := ^^det($,(xj))i<i_^<^. 

The subset of consisting of all finite energy Slater determinants is noted It holds 

that VVf= Wf""®, and, VVf^ for N>2. 

For a wave-function e we define the corresponding A^-body density matrix Ftjj = 

|d')(^'|, which corresponds to the projection on {C^'} in A^i L^(R^,C^). The set of pure-state 
(resp. Slater-state) A^-body density matrices is 

G^"''® := {Ftp, T 6 resp. G^"^*®^ := {Ftp, d- 6 . 
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It holds that = GP'^"® and that ^ G^"™ for iV > 2. The set of mixed-state A^-body 

density matrices G™'^®*^ is defined as the convex hull of G^"^®: 

{ oo oo ^ 

E nk\'i’k)'i>kl 0 < Ufc < 1, E = 1: 6 wr" ■ 

k=l k=l ) 

It is also the convex hull of G^'^*®''. The kernel of an operator F 6 G™'™*^ will be denoted by 

r(ri,si,---,rAr,S7v;ri,Sij---jr^,s^). 


The quantities of interest in density-functional theory are the spin-density 2x2 matrix, and 
the paramagnetic-current. For F € G™’'®'^, the associated spin-density 2x2 matrix is the 2x2 
hermitian function-valued matrix 


i?r(r) := 



(r), 


where, for a,/3 e 

pfi^)-=N E £3,;v 

In the case where F comes from a Slater determinant we get 


N 


Rr{r) = E 

k=l 




^Wk\ 

\<Pi\V' 


( 1 ) 


The total electronic density is pr = Pp + Ppt and the spin angular momentum density is mr = 
trc 2 [cri?r], where 


O' - i^Ox 1 Oy J (7 



contains the Pauli-matrices. Note that the pair (pr,mr) contains the same information as i?r, 
hence the A^-representability problem for the matrix R is the same as the one for the pair (p, m). 
However, as noticed in [7], it is more natural mathematicaly speaking to work with Rr- The 
Slater-state, pure-state and mixed-state sets of spin-density 2x2 matrices are respectively defined 
by 


J Slater . f td t~' ^ /^Slater 1 

N ■= t 6 I ! 

jr^:= {i?r, reG^I, 

J mixed . ( td t~' ^ /^mixedl 

N ■= 1-f^r, t 6 I ■ 


Since the map F i-i- Rr is linear, it holds that c c that is convex, and 

is the convex hull of both and 


For a 7V-body density matrix F 6 GJ^'^'®'^, we define the associated paramagnetic current jr = 
Jr +jr ’^itb 


j“ = Im(7V E f Vr'Pf 


r, a,z,s;r',a,z,s)| dz 


)L. 


In the case where F comes from a Slater determinant iS[$i,...,‘I>Ar], we get 


jr= Elm + 

k=l 


( 2 ) 


Note that while only the total paramagnetic current j appears in the theory of C(S)DFT, the pair 
is sometimes used to design accurate current-density functionals (see [5] for instance). In 
this article however, we will only focus on the representability of j. 
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3 Main results 


3.1 Representability in SDFT 

Our first result concerns the characterization of and For this purpose, we 

introduce 


C^r:={i?6Al2x2(i^(K^C)), = i?>0, 

f trc2[i?] = fV, \/i?6^f2x2(i^HK^C))|, 

Jk3 ) 


(3) 


and := {R € Cn, det R = 0}. The following characterization of Cn was proved in [7]. 


Lemma 1. A function-valued matrix R = 


Ip' cr\ . . 

[a 


in Cm iff its coefficients satisfy 

p^ = N, 


p^!^ > 0 , pV'*' “ kP ^ 0i f ^ f 

a, Vdet(i?) € VF1’3/2 (r3) 

\yafp-^€L\R^), 

|vVdet(i?)|%“^ eL^R^). 


(4) 


The complete answer for 7V-representability in SDFT is given by the following theorem, whose 
proof is given in Section 14.11 


Theorem 1. 

Case N = 1: R holds that 


^Slater ^ ^pure ^ ^0 


Case N > 2; For all N >2, it holds that 


J Slater _ /ypure _ /r-mixed _ p 
N ~ N ~^N- 

Note that the equality for all fV 6 N* was already proven in [7]. 

Remark 1. Gilbert m, Harriman m and Lieb W proved that the N-representability set for 
the total electronic density p is the same for Slater-states, pure-states and mixed-states, and is 
characterized by 

In := {p e p>0, f^^p= N, ■ (5) 

Comparing and we see that our theorem is a natural extension of the previous result. 


3.2 Representability in CSDFT 

We first recall some classical necessary conditions for a pair (R, j) to be fV-representable (we refer 
to [TTl [T2] for the proof). In the sequel, we will denote by p^ := p‘^f p^ := p'*"*' and cr := p^'*' the 

elements of a matrix R, so that R = ^ j, and by p = p^ + p'*' the associated total electronic 

density. 

Lemma 2. If a pair (i?, j) is representable by a mixed-state N-body density matrix, then 


R e Cm 

|jP/p6il(R3). 


( 6 ) 
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From the second condition of m, it must hold that the support of j is contained in the support 
of p. The vector v := is called the velocity field, and w := curl(v) is called the vorticity. 


Let us first consider the pure-state setting. Recall that in the spin-less setting, in the case 
TV = 1, a pair (p,j) representable by a single orbital generally satisfies (provided that the phases 
of the orbital are globally well-defined) the curl-free condition curl(p"^j) = 0 (see [121 lll| l. This 
is no longer the case when spin is considered, as is shown is the following Lemma, whose proof is 
postponed until Section 

Lemma 3 (CSDFT, case N = 1). Let $ = ((/)^, (()■*')^ e yyp^ter such that both (j)^ and (j)^ have well- 
defined global phases in (^^(K.). Then, the associated pair (i?,j) satisfies R e Ci, |jp/p e L^(R^), 
and the two curl-free conditions 


curl 


(i 


Im(crVg) \ ^ ^ 
ppJ- ) 



Im(crVg) \ ^ ^ 
PP^ ) 


(7) 


Remark 2. If we write a = |cr|e''^, then, jcrp = pVS 

Im (ctVct) = IcrpVr = p^X/r. 


( 8 ) 


In particular, it holds that 


curl I 


Im((TV(T) Im((TVcr)\ 


PP'’ 


PP' 


= curl (Vr) = 0, 


so that one of the equalities in o implies the other one. 


Remark 3. We recover the traditional result in the spin-less case, where cr e 0. 

In the case > 1, things are very different. In [T^, the authors gave a rigorous proof for the 
representability of the pair (p,j) by a Slater determinant (of orbitals having well-defined global 
phases) whenever > 4 under a mild condition (see equation @ below). By adapting their proof 
to our case, we are able to ensure representability of a pair (i?,j) by a Slater determinant for 
A^ > 12 under the same mild condition (see Section IT7II1 for the proof). 

Theorem 2 (CSDFT, case N > 12). 

A sufficient set of conditions for a pair (R, j) to be representable by a Slater determinant is 
• Re Cn with N > 12 and j satisfies |jp/p e L^(]R^) 


• there exists <5 > 0 such that, 

sup /(r)^^^'^^^^|w(r)| < oo, sup /(r)^^’^‘^^^^|Vw(r)| < oo, (9) 

reIR3 reIR3 


where w := curl {p ^j) is the vorticity, and 


fir) :=il + ir,)^)il+ir,)^)il^ir,)^). 


Remark 4. The eonditions are the ones found in The authors conjectured that this 

condition "can be considerably loosened". 

Remark 5. IFe were only able to prove this theorem for N > 12. In JT^, the authors proved that 
eonditions were not sufficient for N = 2. We do not know whether conditions ® are sufficient 
in the case 3 < A'^ < 11. 

Let us finally turn to the mixed-state case. We notice that if (i?, j) is representable by a Slater 
determinant J^[<I>i,..., 4*^?], then, for all k e N*, the pair (A:/A^)(R, j) is mixed-state representable, 
where N is the number of orbitals (simply take the uniform convex combination of the pairs 
represented by c5^[$i], .5^[$2], etc.). In particular, from Theorem [21 we deduce the following 
corollary. 
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Corollary 1 (CSDFT, case mixed-state). 

A sufficient set of conditions for a pair {R,j) to be mixed-state representable is R € C% for some 
N 6 N*, j satisfies |jp/p e and ([5]) holds for some (5 > 0. 

In the authors provide different sufficient conditions than ([2]) for a pair (p, j) to be mixed- 
state representable, where p is the electronic density. They proved that if 

(l + |-p)p|v(p-'j)|%L'(K3), 

then the pair (p, j) is mixed-state representable. Their proof can be straightforwardly adapted for 
the representability of the pair so that similar results hold. The details are omitted here 

for the sake of brevity. 


4 Proofs 

4.1 Proof of Theorem [1] 

The mixed-state case was already proved in [7] . We focus on the pure-state representability. 


Case = 1 

The fact that simply conies from the fact that = GP""®. To prove jfiater ^ 

C°, we let R e jySiater represented by $ = e so that 

W \^fl' 

Since R e jySiater ^ j^-mixed ^ R=0,we deduce ReCi. 

We now prove that c jySiater j^g|. ^ ^ ^ j g deti? s 0 and Lemma [H we get 

fp^/^>o, pV = kP, f [ P^ = i, 

JR3 Jr3 

^Vaf/p€L\R^). 

There are two natural choices that we would like to make for a representing orbital, namely 



Unfortunately, it is not guaranteed that these orbitals are indeed in H^{R^,C^). It is the case 
only if IVcrp/p'*’ is in for <i)i, and if |Vcrp/p^ is in ffi(R^) for $ 2 - Due to (fT0)l . we know 

that |Vcrp/p € L^(]R^). The idea is therefore to interpolate between these two orbitals, taking $i 
in regions where p^ » pp and <I >2 in regions where p^ » p^ This is done via the following process. 


Let X s G°“(]R) be a non-decreasing function such that 0 < y < 1, x(^) 
xix) = 1 if a: > 1. We write cr = a + i/3 where a is the real-part of a, and /3 is 
We introduce 


Ai 

A 2 


^JoP^TxRjf[pi)ffi 

ciAi + /3pi 


Pi 

P2 


■■= Vi-x^(pVp') 

l3Xi-api 



= 0 if a; < 1/2 and 
its imaginary part. 


and we set 


At - 


:= Ai + ipi and /)■*' := A 2 + ip 2 - 


6 











Let us prove that $ represents R and that $ := ^ First, an easy calculation shows 

that 

,.t ,2 ^ 2 , 2 + kP t 

W\ =Ai+/ri =- 

|,P2 + ^ i 

Re = A 1 A 2 - P 1 P 2 = p = a. 


Im 


P' 

/3(A^ + /^f) 


- A1P2 + A2P.1 - 


so that $ 6 L^(K^,C^) with ||$|| = 1, and $ represents R. To prove that $ e ypp^ter^ need to 
check that Ai,A 2 ,pi and ^2 are in iL^(R^). For Ai, we choose another non-increasing function 
^ 6 (^“(R) such that 0 < ^ < 1, ^(x) = 0 for x < 1, and ^(x) = 1 for x > 2. Note that (1 - x)^ = 0. 
It holds that 

VAi = (1 - f (pVk)) VAi + e"(pVk) VAi. (12) 

The second term in the right-hand side of (HID is non-null only if > p^, so that xipVP^) = 1 on 
this part. In particular, from the equality p^p^ = |tTp, we get 

^^(pVk)Ai= eipVp^)^=eiP^iP^)V?, 

and similarly, 

f (pVp^)VAi = f(pVk)vV7, 

which is in L^(]R^) according to dTIlD . On the other hand, the first term in the right-hand side 
of (fT^ is non-null only if p^ < 2p^, so that (l/3)p < p^ on this part. In particular, from the 
following point-wise estimate 

lvV7l + IvVpI, 

which is valid almost everywhere whenever /,p > 0, the inequality (a -n b)^ < 2{a? + 6^), and the 
fact that + x^/3^ < kP, we get on this part (we write x for x(pVp’*’)) 

2 


|VAip = 


+ y2^2 _ y/gj2 + ^2^2 


\ p^ {p^r 


< 2 


p^Vpj-p’Vp* I 


|VQ!p_^^kA (pl )2 


, 2x^|V/3p ^ 2H 


(k) 


HvVkpj- 


We hiially use the inequality {p^) ^ < {2>lp), and the inequality kP/(pk^ = pVp^ - 2 and get 

^ |r7^,i2 /Ir7^t|2 lr7^J-|2’' 


ivAy<c(i™ 


|V/3p 


ivvFfV 


The right-hand side is in L^(K^) according to (fTUl) . Hence, (1 -‘?^(pVpk) |VAi| 6 L^(R^), and 
finally 

The other cases are treated similarly, observing that. 
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• whenever , then x = 1 ) $ = $i where $i was defined in (ED- We then control 

with the inequality {p‘^)~^ <2p~^ ; 

• whenever p'' < p^ 12, then y = 0, $ = $ 2 - We control {p^)~^ with the inequality {p^)~^ < |p”^ 


whenever p ^/2 < p^ < p\ then both (p^) ^ 

{p^)-^<2p-\ 


and p^ are controlled via (p^) ^ < 3p ^ and 


Case N >2. 

Since J'®'’"*®"' c c = C^v, it is enough to prove that Cn c We start with a key 

lemma. 


Lemma 4. For all M,N it holds that jf+M = 

Proof of Lemma\^ The case is trivial: if i? € is represented by 

the Slater determinant oS^[<i)i,... ^at+m], then, by denoting by i?i (resp. R^) the spin-density 
matrix associated to the Slater determinant at] (resp. c5^[$Ar+i,..., d)Ar+i]), it holds 

R = Ri + R 2 (see Equation ([T]) for instance), with i?i 6 and R 2 € 


The converse is more involving, and requires an orthogonalization step. Let i?i e 
be represented by the Slater determinant and R 2 e be represented by 

the Slater determinant o5^[$i,...,$m]- We cannot directly consider the Slater determinant 
^[$i,...,$Ar,$i,...,$M], for ($i,...,$ 7 v) is not orthogonal to ($i,...,$m)- 


We recall the following lemma, which is a smooth version of the Hobby-Rice theorem m (see 
also uni), and that was proved by Lazarev and Lieb in [T3] (see also IE)- 

Lemma 5. For all N e N*, and for all (/i,...,/a/) e L^(IR^,C), there exists a function u e 
(^“(K^), with bounded derivatives, such that 

Vl<fc<iV, [ /fee‘“ = 0. 

Jb3 

Moreover, u can he chosen to vary in the ri direction only. 

We now modify the phases of $ 1 ,..., as follows. First, we choose ui as in Lemma [5] such 
that, 

Vl<fc<iV, f^^{7Jl+7A)e^ = 0, 

and we set $a/+i = $ie'“b Note that, by construction, $jv+i is normalized, in iJ^(IR.^,C^), and 
orthogonal to ($i,...,$jv)- We then construct U 2 as in Lemma [S] such that 

yi<k<N + l, +44)e‘'^ = 0, 

and we set ^n +2 = We continue this process for 3 < fc < M and construct $A/+fe = 

We thus obtain an orthonormal family ($ 1 ,..., ^at+m)- By noticing that the spin-density matrix 
of the Slater determinant is the same as the one of ■ ■ ■ ,^n+m] (the 

phases cancel out), we obtain that R = Ri + R 2 , where R is the spin-density matrix represented 
by .., $Af+M]- The result follows. 

□ 

We now prove that Cat c JT^'®*®''. We start with the case N = 2. 


Case N = 2. 

Let R = 6 C 2 . We write \/R = j, with r\r^ e (i7^(]R^))^ and s in 7J^(K^,C). Let 


It is easy to check R = + R^, that R^^^ are hermitian, of null determinant, and 6 

■Ad 2 x 2 C)). However, it may hold that /R 3 trc 2 [i?^] ^ N*, so that R^ is not in for 

some M € N*. 


The case = 0 or !?■*■ = 0 are trivial. Let us suppose that, for a e {t, J.}, m“ := pRc + 0. In 
this case, the matrices are in C°, hence are representable by a single orbital, due to 

the first statement of Theorem[T] Let $ = e and $2 = (</> 2 i'^ 2 )^ ^ 

be normalized orbitals that represent respectively R^ and It holds 

= {m^y^R^ and $ 2^2 = ^ = 


From the Lazarev-Lieb orthogonalization process (see Lemma[5]), there exists a function u € C'°“(R) 
with bounded derivatives such that 

($i|$2e'“> = e™ = 0. (14) 


Once this function is chosen, there exists a function v € C°“(]R) with bounded derivatives such 
that 


= ($i|$2e‘(“^’')) = ($2e™|$ie‘") = ($2|$2e'") = 0. (15) 


We finally set 



From (fH)) . we deduce ||4>ip = ||<I> 2 p = 1, so that both $1 and ^2 are normalized. Also, from dT51) . 
we get (4 )i|<I> 2) = 0, hence {$ 1 , 4 ) 2 } is orthonormal. As $1 and $2 are in iL^(]R^,C^), and u and 
V have bounded derivatives, $1 and $2 are in iL^(R^,C^). Finally, it holds that 


$ 1 ${ +$ 2$2 = 

= i + m^$2$2 + 2V rri^m^Yle 

+ + m'*'$2$2 “ 2'/mJmJRe ($i$26~™) j 

= + TO^$2$2 = R- 


We deduce that the Slater determinant =5^[$i,$2] represents R, so that R 6 Altogether, 

C 2 C jSiater^ therefore C 2 = 


Case N > 2. 

We proceed by induction. Let R e Cn+i with N > 2, and suppose Cn = We use again 

the decomposition m and write R = R^ + R^, where R^^^ are two null-determinant hermitian 
matrices. For a € {t, I}, we note m“ := pRa. Since + 1 > 3, at least m} or m} is 

greater than 1. Let us suppose without loss of generality that > 1. We then write R = Ri + R 2 
with 

i?i := (m^)~^i?^ and R 2 ■= R^ + R^) . 

It holds that i?i 6 C® = j/'^iater R 2 e Cn = induction). Together with Lemma SJ we 

deduce that R € The result follows. 
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4.2 Proof of Lemma [3] 

Let $ = having well-defined global phases in C^(M), and let be the 

associated spin-density matrix and paramagnetic current. It holds 



For a € {t, !•}) we let t“ be the phase of (/)“, so that (j)°‘ = . Setting r = r'^ - r^, we obtain 

cr = |cr|e*’^ = The paramagnetic current is 

j = -I- p'^'Vt'*' = pVr’*' -I- p'^ Vt = pVr^ - p'*'V'r. 


In particular, using ([5]), 

j Im (ctVct) j-p^Vr 

p ppJ- p 

is curl-free, and so is 

j ^ Im(aVg) ^ 

p ppt 


4.3 Proof of Theorem [2] 

We break the proof in several steps. 


(16) 


(17) 


Step 1: Any R e Cn can be written as R = Ri + R 2 + R 3 with Rk 6 Nk > 4. 

Let R = 6 Cat, with N > 12. We write \/i? = with e and s in 

i7^(R^,C). We write R = where R^^^ were defined in (fT^ . As in the proof of Theorem [1] for 

the case N = 2, R}^^ are hermitian, of null determinant, and c A42x2 (i7^(K.^, C)). However, 

it may hold that f trc 2 [i?'’^] N*, so that R^ is not in for some M eN*. In order to handle 

this difficulty, we will distribute the mass of R^ and R^ into three density-matrices. 

More specifically, let us suppose without loss of generality that f trc 2 [i?'^] > f trc 2 [i?'^]. We set 


i?i = (i-ei)i?^ + 6-R^ 

i?2 = 6(l-6)i?^ (18) 

i?3 = (l-6)l?^ + 6-R^ 


where are suitable non-decreasing functions in (^“"(R^), that depends only on (say) ri, 

such that, for 1 < fc < 3, 0 < < 1. We will choose them of the form ^fe(r) = 0 for Xi < ak and 

^fc(r) = 1 for all xi > I3k > ctk, and such that 

(1 - ^ 1)6 = (1 - 6)^3 = (1 - Ci)C3 = 0. (19) 

Finally, these functions are tuned so that f trc^{Rk) € N* and f tic'^{Rk) > 4 for all 1 < fc < 3 
(see Figure [T] for an example of such a triplet (^ 1 ,^ 2 , ^ 3 ))- Although it is not difficult to convince 
oneself that such functions ^k exist, we provide a full proof of this fact in the Appendix. 

From 119]), it holds that for all 1 < fc < 3, i?fc e , and that Ri + R 2 + R 3 = R^ + R^ = R. 

Step 2 : The pair (i?i,ji) is representable by a Slater determinant. 

In order to simplify the notation, we introduce the total densities of R^ and i?b 

/^ := |r^P-t-|sp and /^ := + |sp. 
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Figure 1: Weights of the matrices (black) and (gray) in (a) i?i = (1 - ^i)R^ + ^ 2 RK (b) 
i?2 = 6(1 and i?3 = (1 -6)-R^ + 61?^- 


Recall that p = + f^. We consider the previous decomposition R = Ri + R 2 + R 3 , and we 

decompose j in a similar fashion. More specifically, we write j = ji + j2 + js with 


ji = (i-6) 


(^j - Im (Ws) j + 6 + Im (Ws) j , 

h = 6(1-6) -Iin(Ws)j , 

ja = (1-6) +lm(sVs)j +6|^j -lm(Ws)j . 


( 20 ) 


Let us show that the pair (i?i,ji) is representable. Following [T^], we introduce 


1 

6^) = - / 
m J- 


rdy, 


/-OC (l + y2)(l+5)/2 

where <5 is the one in ([HI), and m is a constant chosen such that 6°°) = 1- We then introduce 

2 

?7i,i(r) = —6r + a), 


?7i,2(r) = + ^)(1 - ?7i(r)), 

2 

?7i.3(i-) = -^^^ 62^2 +7)(l-?7i(r) -’72(r)), 


( 21 ) 


7i.fc(r) = 


1 


-(1-? 7 i(r)-? 72 (r)- 773(r)) for 4 < fc < iV, 


N-3 

where a,/3,7 are tuned so that, if pi := trc2i?i denotes the total density of Ri, 

VI < /c < Nk, f m,kPi = 1- 

JR3 

It can be checked (see [12]) that 771 fc > 0 and T,k=i 7i,fc = 1- We seek orbitals of the form 

:= j + V 6 (j) 1 < fc < w, 


(22) 
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and where the phases ui^k will be chosen carefully later. From m, we recall that (1 - ^ 1)^2 = 0, 
so that, by construction, $1 ^ is normalized, and 


= ViykRi- 

Let us suppose for now that the phases ui^k are chosen so that the orbitals are orthogonal. This will 
indeed be achieved thanks to the Lazarev-Lieb orthogonalization process (see Lemma [51). Then, 
7 v] indeed represents the spin-density matrix Ri. The paramagnetic current 
of 4* is (we recall that and are real-valued, and we write s = for simplicity) 

Ni 

- Cl) {\r^fyui,k + |sp V(-T + Upfc)) + 


k=l 


Ni 


+ Z Vi,k^2 (|sp v(r + + |r^p Vui,fc) 


k=l 


= ((1 - Ci)/^ + C 2 /O ( E Vui,fe j + (6 - (1 - Cl)) kp Vr. 

Since |spVT = Im(sVs), this current is equal to the target current ji defined in (EUl) if and only if 


Ni 


Pi- = piYj i?fcVui,fc- 
p 


(23) 


k=l 


In [12], Lieb and Schrader provided an explicit solution of this system when iVi > 4. We do not 
repeat the proof, but emphasize on the fact that because condition ([5]) is satisfied by hypothesis, 
the phases ui^k can be chosen to be functions of ri only, and to have bounded derivatives. In 
particular, the functions are in iJ^(R^,C^). Also, as their proof relies on the Lazarev-Lieb 
orthogonalization process, it is possible to choose the phases ui^k so that the functions $1 ^ are 
orthogonal, and orthogonal to a finite-dimensional subspace of L^(]R^,C^). 


Altogether, we proved that the pair (i?i, ji) is representable by the Slater determinant ^[$ 11 ,..., $1 jvi]- 


Step 3: Representability of (i? 2 j 2 ) and (i? 3 ,j 3 ), and finally of (i?,j). 

In order to represent the pair (i? 2 ,j 2 )) we first construct the functions 772 ,fe for 1 < fc < AI 2 of the 
form (l?I|) so that holds for p 2 '■= trc 2 i ?2 • We then seek orbitals of the form 

d’ 2 ,fe := for l<k<N 2 . 

Reasoning as above, the Slater determinant of these orbitals represents the pair (R 2 J 2 ) if and 
only if 

h 

P2— = P2 E 172,fcVM2.fc- 
p k=l 

Again, due to the fact that N 2 > 4, this equation admits a solution. Moreover, it is possible to 
choose the phases U 2 ,k so that the functions $ 2 ,fe are orthogonal to the previously constructed d*! fc. 

We repeat again this argument for the pair (i? 3 ,j 3 ). Once the new set of functions 773 .^ is 
constructed, we seek orbitals of the form 

^>3.;= := -6) j + j j 

and construct the phases so that the functions 4)3 fc are orthogonal to the functions d*! fc and ^ 2 ,k- 

Altogether, the pair (i?, j) is represented by the (finite energy) Slater determinant o5^[d)i 1 ,..., d^i^Ar^, d) 2 ,i,..., $ 2 ,at 
which concludes the proof. 
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5 Appendix 

We explain in this section how to construct three functions s like in Figure [TJ 

In order to simplify the notation, we introduce 

/(^) ■= ff trc2(i?'‘)(r,r2,r3) dradra, 

JjRxR 

g{r)-= II trc 2 (i?^)(r,j/, 2 :) dradrg, 

JjRxR 

where were defined in (HU. We denote by 

F{a) = f f(x)dx and G(a) = f g{x)dx, 

J—OO sj—oo 

and finally T = F{oo) = f^f and G = G(oo) = f^^g. Note that F and G are continuous non¬ 
decreasing functions going from 0 to (respectively Q), and that it holds T + G = N. Let us 
suppose without loss of generality that F < G, so that 0<F<Nj2<G<N. If.7^ = 0 , then 
!?■*■ = 0 and we can choose Ri = R 2 = {dlN)R} e C 4 and R^ = {N - 8)INR^ 6 C^_g. Since N > 12, 
it holds A ^-8 > 4, so that this is the desired decomposition. We now consider the case where F + 0. 

In order to keep the notation simple, we will only study the case .F < 8 (the case .F > 8 is 
similar by replacing the integer 4 by a greater integer M such that F < 2M < - 4 in the sequel). 

We seek for a such that 

I iZo fix)dx < 4 and f{x) + g{x) > 4, 

1 /a”< 4 and fZ^g{x)dx + f^f{x)dx>4:, 

or equivalently 

J F(a) < 4 and F(a) + G - G(a) > 4, 

\F-F{a)<d and F - F{a) + G{a) > F 

that is 

JP-4<F(a)<4, and F{a) + d-F < G{a) < F{a) + G - d. (24) 

Let be such that E'(a(j;r_ 4 )) = F -d (with a(j;r_ 4 ) = -00 if .F < 4), and a( 4 ) be such 

that F(a( 4 )) = 4 (with 0 ( 4 ) = +00 if < 4). As F is continuous non-decreasing, the first equation 
of (IMl) is satished whenever a(jF- 4 ) < a< a( 4 ). 

The function [Q!(^_ 4 ),a 4 ] 9 a 1 -^ m(a) := F{a) +d-F goes continuously and non-decreasingly 
from 0 to 8 - .F, and the function ^ a M{a) := F{a) + G - d goes continuously 

and non-decreasingly from N - 8 to G between a(^_ 4 ) and 0 ( 4 ). In particular, since G{a) goes 
continuously and non-decreasingly from 0 to G, only three cases may happen: 


• There exists ag s ('a(:r- 4 ),Q!( 4 )) such that m(ao) < G(ao) < M(ao). In this case, (1^ holds 
for a = ao. By continuity, there exists e > 0 such that 

f F{a + e) <d, 

-I F(a)-I-^ - G(q! + e) > 4, 

[ G(a) -HF'-F(a + e) >4. 
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Let ^2 e (^“(R) be a non-decreasing function such that ^ 2 ( 2 ;) = 0 for a; < a and ^ 2 ( 2 ^) = 1 for 
x> a + £. Then, as 0 < ^2 ^ 1, it holds that: 


and 


f < F{a + £) <4 

JR 

Ai-6)/+ /" g> F{a)+g-G{a + £)> 4 . 
JR J a+s 


such that ^ 3 (x) = 0 for 


such that ^ 1 ( 2 ;) = 1 for 


We deduce that there exists an non-decreasing function ^3 e 
X < a + £, and such that 

/"(I - 6)/ + ^39 = 4. 

Jr 

Note that (1 -^ 2)^3 = 0. On the other hand, from 
f 46/<^-n«)<4 
1 fR^2f + > F - F{a +e) + G{a) > 4, 

we deduce that there exists an non-decreasing function e 
x> a, 

[ (1-Ci)5 + ?2/ = 4. 

Jr 

and (1 -^ 1)^2 = (1 -^ 1)^3 = 0- Finally, we set 

Ri = il-^i)RU^ 2 R^ 

R 2 = ^i{l-i3)R' 

By construction, R = R'^ + R^ = Ri+ R 2 + R^, Ri e and R 3 eC^. We deduce that R 4 e . 
Together with the fact that N > 12, this leads to the desire decomposition. 

For all a e (a(:r- 4 ), q;( 4 )), it holds G(a) < m(a). Note that this may only happen if TO(a( 4 )) > 
0, or < 4, so that ^ > 7V-4 > 8. It holds G(a(;r- 4 )) = 0, so that g(r) is null for r < Q!(^_ 4 ). 
Let ao be such that a(;r- 4 ) < ao < ck( 4 )- As 


f f = R>4 and f"' = T-F{ao)<4, 

JR J Q;o 


there exists a non-decreasing function e G°°(R) satisfying ^i(x) = 1 for a; > ao and such 
that 

/ 6 / = 4. 

Jr 

Now, since G(a( 4 )) < m(a( 4 )) = 8 - J^, it holds that 
f /R(l-ei)/<A(a(4))=4 

1 /fl(l - 6 )/ + > A(a(^_4)) + 0 - G(a(4)) > 4. 

There exists a non-decreasing function ^2 e G““(M) satisfying ^ 2 ( 2 ;) = 0 for x < ao and such 
that 

f (1 “Cl)/ + C 2 J = 4. 

Jr 

Note that (1 -^ 1)^2 = 0. Finally, we set 

Ri = ^iR^ 

i? 2 =(l-C 2 )i?^ 

i?3 = C 2 i?^ + (l-Cl)^^- 

By construction, it holds that R = Ri + R 2 + R^, and that i?i e C 4 and R 3 e C 4 . We deduce 
R 2 e Cn- 8 ’ result follows. 
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• For all a € (q:(^_ 4 ),< 2 ( 4 )), it holds G(a) > M(a). This case is similar than the previous one. 
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